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ABSTRACT 



We present near-infrared (1 - 2.5 ^m) spectroscopic and photometric results of 
Nova V2615 Ophiuchi which was discovered in outburst in 2007 March. Our observa- 
tions span a period of ~ 80 days starting from 2007 March 28 when the nova was at 
its maximum Hght. The evolution of the spectra are shown from the initial P-Cygni 
phase to an emission-line phase and finally to a dust formation stage. The charac- 
teristics of the JHK spectra are very similar to those observed in a nova outburst 
occurring on a carbon-oxygen white dwarf. We analyse an observed line at 2.088 /zm 
and suggest it could be due to Fell excited by Lyman a fluorescence. The highlight of 
the observations is the detection of the first overtone bands of carbon monoxide (CO) 
in the 2.29 - 2.40 /im region. The CO bands are modeled to estimate the temperature 
and mass of the emitting CO gas and also to place limits on the ^^C/^'^C ratio. The 
CO bands are recorded over several epochs thereby allowing a rare opportunity to 
study its evolution from a phase of constant strength through a stage when the CO is 
destroyed fairly rapidly. We compare the observed timescales involved in the evolution 
of the CO emission and find a good agreement with model predictions that investigate 
the chemistry in a nova outflow during the early stages. 

Key words: infrared: spectra - line : identification - stars : novae, cataclysmic vari- 
ables - stars : individual (V2615 Oph) - techniques : spectroscopic 



1 INTRODUCTION 

V2615 Ophiuchi (Nova Ophiuchi 2007) was discovered by 
Nishimura (2007) at a visual magnitude of 10.2 on 2007 
March 19.812 UT. No star was visible at the position of the 
nova on a film taken by Nishimura on 2007 March 17.82 
UT and earlier survey films since 2005 (limiting mag 11.5- 
12). The discovery report was supplemented by several other 
observers from observations obtained close to and around 
the time of discovery. These reports, in conjunction with the 
lack of a detection of the nova in images taken just a few 
days prior to discovery, imply that the nova was discovered 
a few days before maximum light which was reached on 28 
March. 

The early low-resolution optical spectra (410-670 nm, 
resolution 1500) of V2615 Oph by Naito and Narusawa 
(2007) on March 20.84 UT showed Balmer-hnes having P- 
Cyg features and Fe II lines (multiplets 42, 49, 74), suggest- 
ing that the nova is a "Fe H"- type nova. The FWHMs of the 
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Hq, H/3, and H7 emissions were 920, 810, and 790 km/s, and 
the displacement of the P-Cyg absorptions from the Balmer 
emission peaks were 940, 820, and 830 km/s, respectively. 
Early optical spectra were also obtained by Munari et al. 
(2007, 2008) on Mar. 22.17 and 24.18 UT. On March 22.17 
UT, the spectrum was characterized by weak emission lines 
with P-Cygni profiles. The strongest emission lines were due 
to Fe II multiplets 27, 28, 37, 38, 42, 48, 49, 55, and 74; Si II 
multiplets 2, 3, 4 and 5; N I multiplet 3; Ca II H and K; Na 
I Dl and D2; and O I 777.2-nm (the O I line showed the sec- 
ond strongest emission component after Ha). With respect 
to the peak of the Ha emission, the absorption terminal ve- 
locity was -1415 km/s, and the mean velocity -910km/s. In 
comparison to the March 22.17 UT spectrum, the Mar. 24.18 
UT spectra show large changes: mainly a marked reduction 
in equivalent width of emission lines and reduction in out- 
flow velocity of absorption components of P-Cyg profiles. 

Infrared observations of V2615 Oph have been made by 
Das et al. (2007) and Rudy et al.(2007) - some of the findings 
in Rudy et al. (2007) are discussed in details subsequently. 
The highlight of the Das et al. (2007) report was the de- 
tection of first-overtone CO emission in the nova - similar 
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emission has been recorded previously in only a few novae. 
In this paper, we present JHK spectrophotometry of V2615 
Oph between March 28 and June 9. The spectroscopic re- 
sults are fairly extensive as spectra have been obtained on 
15 occasions during this period. 



2 OBSERVATIONS 

The Near-Infrared JHK spoctrophotometric data of V2615 
Oph presented here were acquired at the Mount Abu 1.2m 
telescope. The first observation was recorded on the night of 
2007 March 28 while the nova was close to the optical maxi- 
mum. The nova was then observed regularly until the arrival 
of the monsoon, when inclement weather brought our moni- 
toring to an end. The instrument used was the Near Infrared 
Imager/Spectrometer with a 256x256 HgCdTe NICMOS3 
array. The observation logs for spectroscopy and photom- 
etry are presented in Tables 1 and 2 respectively. Near-IR 
JHK spectra were obtained at similar dispersions of ~ 9.75 
A/pixel in each of the J, H, K bands. In each band, spectra 
were recorded at different positions along the slit (slit width 
~1 arcsec) to provide for the background subtraction. To 
remove telluric features in the spectra of V2615 Oph, we 
obtained spectra of a nearby standard star (SAO 184301; 
spectral type AOV) at a similar airmass as the target. 

The spectra were extracted using the APEXTRACT 
task in IRAF and the spectra wore calibrated in wavelength 
using a combination of OH sky lines and telluric lines in 
the extracted spectra. Following the standard reduction pro- 
cedure, the nova spectra were then ratioed with the spec- 
tra of the standard star from which the hydrogen Paschen 
and Brackett absorption lines had been manually removed. 
These ratioed spectra were multiplied by a blackbody curve 
corresponding to the effective temperature of the standard 
star to yield the final spectra. 

V2615 Oph was monitored photometrically in the JHK 
bands for 85 days after discovery. Photometric observa- 
tions were performed under photometric sky conditions us- 
ing the imaging mode of the NICMOS3 array. In each of 
the J, H, K filters, several frames at 5 dithered positions 
offset typically by 20 arcsec, were obtained of both the 
nova and a selected standard star (2MASS J16232693 - 
2425291; J=7.340,JJ=6.027, 7^=5.464). Near-IR JHK mag- 
nitudes were then derived using the IRAF aperture pho- 
tometry task APPHOT. The derived JHK magnitudes are 
given in Table 2, the typical errors associated with them lie 
in the range of 0.02 to 0.04 magnitudes. 



3 RESULTS 

3.1 Optical and nesir-infrared lightcurve 

The optical lightcurve of V2615 Oph is presented in Figure 
1. The object was discovered about 9 days before its max- 
imum which was reached on 2007 March 28 at V = 8.52. 
On the whole, the object shows a steady post-maximum de- 
cline in brightness of 0.05 mag/day for the first 80 days 
after maximum. Subsequently there is a steep decline in the 
hghtcurve due to dust formation - this aspect is addressed 
shortly while discussing the near-IR light curve. From the 
lightcurve we estimate t2 and ts - the time for a drop of 2 



Table 1. A log of the spectroscopic observations of V2615 Oph. 
The date of outburst is taken to be 2007 March 19.812 UT 



Date Days Integration time 

2007 since (sec) 

(UT) Outburst JHK 



Table 2. A log of the photometric observations of V2615 Oph. 
The date of outburst is taken to be 2007 March 19.812 UT 



Date 


Days 




Magnitudes 




2007 


since 








(UT) 


Outburst 


J 


H 


K 


Mar. 28.993 


09.181 


6.26 


5.83 


5.39 


Mar. 31.961 


12.149 


6.94 


6.51 


5.90 


Apr. 02.970 


14.158 


6.49 


6.07 


5.54 


Apr. 05.993 


17.181 


6.54 


6.14 


5.64 


Apr. 07.956 


19.144 


6.93 


6.57 


6.18 


Apr. 15.979 


27.167 


7.36 


7.07 


6.53 


Apr. 16.981 


28.169 


7.22 


6.97 


6.46 


Apr. 18.940 


30.128 


7.30 


6.96 


6.62 


Apr. 26.885 


38.073 


7.76 


7.53 


7.09 


Apr. 30.906 


42.094 


7.98 


7.68 


7.18 


May. 02.932 


44.120 


8.14 


7.84 


7.07 


May. 06.934 


48.122 


8.13 


7.72 


6.79 


May. 08.882 


50.070 


7.94 


7.32 


6.26 


Jun. 09.887 


82.075 


9.31 


7.62 


5.34 



and 3 magnitudes respectively in the visual brightness - to 
be 33 and 58 days respectively. The use of these values of 
t2 and ts in various MMRD relationships (della Valle et al. 
(1990), CapaccioU et al. (1989) Cohen (1988) for t2; Schmidt 
(1957) for tz) leads to closely similar values for the absolute 
visual magnitude - a mean value of Mv = 7.16 ± 0.12 is ob- 
tained. The extinction toward the object can bo estimated 
in two ways. Rudy et al. (2007), using the strength of the 
01 lines, have determined the reddening prior to dust for- 
mation to be E(B-V) = 1.0 (or Av = 3.1 magnitudes). This 
leads to an distance estimate of D = 3.25 kpc to the ob- 
ject. To show that this value of Av is reasonable, we use the 
extinction data of Marshall et al. (2006) which shows that 
Aks is constant at = 0.34-0.36 beyond 3 and upto 12 kpc 
in the direction of V2615 Oph. Thus, assurrung Av/Aks ~ 
11 (Koornneef 1983), a maximum value of Av = 3.85 is sug- 
gested toward V2615 Oph thereby leading to a lower limit 
on the distance oi D = 2.3 kpc. For the present work we 
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Figure 1. The V band lightcurve is siiown at tlie bottom (plio- 
toelectric/CCD measurements shown by filled circles ; visual esti- 
mates by open circles) based on data from lAU circulars, AAVSO 
(American Association of Variable Star Observers) and AFOEV 
databases (Association Prancaise des Observateurs d'Etoiles Vari- 
ables, France). At the top are shown the near-IR J, H and K band 
light curves based on the present observations (J : filled circles, 
H : gray triangles - up; K : black triangles - down.) 



adopt the distance estimate derived from the Rudy et al 
(2007) findings viz. D = 3.2 kpc. 

Our IR observations, begun close to optical maximum, 
show a steady decline for the first 42 days. After this each 
of the J,H,K magnitudes show an increasing trend possibly 
indicating the onset of dust formation. However, no corre- 
sponding sharp decline in the optical hght curve is seen at 
this stage implying that any dust, that may have formed, is 
likely to be optically thin. Rudy et al. (2007) have reported 
0.4 to 2.5 fim spectroscopy of V2615 Oph on May 7 2007 
and 0.8 to 5.5 fim spectroscopy on 2007 May 31. They find 
that the nova evolved dramatically between the two observa- 
tions (i.e between May 7 and May 31) due to the formation 
of dust. The reddening derived from the 01 lines increased 
from E(B-V) = 1.0 to 1.3 in their data. Our last observations 
on 9 June 2007 show a significant infrared excess due to 
emission by the newly formed dust. 

3.2 General characteristics of the JHK spectra 

The lines in the earliest spectrum of 28 March, when the 
nova was at maximum hght, show P-Cygni structure with 
a prominent absorption component. By following the Pa /3 
line for example, it is seen that the P-Cygni absorption com- 
ponent persists for at least 10 days (upto 7 April) indicating 
that mass loss extends over a prolonged period. The mini- 
mum of the absorption component is found to be displaced 
from the emission peak by ~ 950 km/s on 28 March but this 
value changes closer to ~ 1400-1500 km/s for 31 March and 
other subsequent dates upto 7 April. Variations, of a similar 
nature, in the outfiow velocities were also noted by Munari 
et al (2007) on 22 and 24 March - when the nova was on its 
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Figure 2. The Paschen /3 1.2818 fim line is presented on different 
days to show the evolution of the profile. The two dotted lines, 
at left of the line center, are positioned at -950 and -1500 km/s 
respectively to indicate the change in position of the P-Cygni 
absorption minimum with respect to the emission peak (further 
details are in the text of Section 3.2) 



rise to maximum light. Clearly, the mass loss process and 
its kinematics during the early stages is not uniform. 

Mosaics of the the J, H and K band spectra are pre- 
sented in Figures 2, 3, and 4 respectively. The early spectrum 
of V2615 Oph is typical of a carbon-oxygen nova (CO nova) - 
examples of which are V2274 Cyg (Rudy et al. 2003) , V1419 
Aql (Lynch et al. 1995) and V1280 Sco (Das et al. 2008). In 
addition to lines of H, He, N and O, these novae show strong 
lines of carbon e.g in the J band in the 1.16 to 1.18 nm wave- 
length region. In contrast, such carbon lines are weak in the 
spectra of ONeMg novae e.g. in the novae V597 Pup and 
V2491 Cyg (Naik et al. 2009). The IR-based classification of 
V2615 Oph is consistent with its optical classification as a 
Fell type nova (Naito & Narusawa, 2007). Fe II novae are be- 
lieved to be associated with explosions on CO white dwarfs 
(Williams 1992). We have presented the line identification 
in Table 3 but do not show them marked on the spectra 
as a separate figure. However, since the lines are very simi- 
lar to those seen in V1280 Sco, the line-identification figure 
for V1280 Sco could be referred to which contains greater 
details (Das et al. 2008). The most prominent lines in the 
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JHK spectra of V2615 Oph are the Paschen and Brackett 
hydrogen recombination hnes; the Lyman fi and continuum 
fluoresced 01 lines at 1.1287 and 1.3164 ytim; the Hel hnes at 
1.0830 jxm and 2.0585. Strong carbon hnes are seen in the J 
band also in the H band rcdwards of Br 11. Among subtle 
features, the presence of a C I line at 1.6 jiia, which could 
be mistaJien as just another member of the H band Brackett 
series lines, should not be missed. The region between 1.2 to 
1.275 /Ltm contains the complex blend of a large number of 
NI and C I lines that arc often seen in the early spectra of 
CO novae. The presence of several lines of Nal and Mgl ( the 
prominent ones being at 1.1404, 1.5040, 1.5749, 2.2056 and 
2.2084 /im ) arc worth noting as it may be inferred, from 
their presence, that dust will form in a nova. During the 
analysis of V1280 Sco, one of the conclusions that emerged, 
was that the lines of Nal particularly, and also of Mgl, are as- 
sociated with low excitation and ionization conditions (Das 
et al. 2008). Such conditions necessarily imply the existence 
of a cool zone which is conducive for dust formation. This 
was observationally corroborated in the case of several no- 
vae, in which these lines were detected, and which proceeded 
to form dust (e.g. V2274 Cyg (Rudy et al. 2003), V1419 Aql 
(Lynch, ct al. 1995), NQ Vul (Ferland et al 1979), V705 
Gas (Evans et al 1996), V842 Cen (Wichmann et al. 1991), 
V1280 Sco (Das et al. 2008); for a more detailed discussion 
sec Das et al. 2008). The presence of these lines in V2615 
Oph, and the dust formation witnessed subsequently, is con- 
sistent with this scenario. The formation of dust is clearly 
reflected in the last spectrum of 9 June. Each of the spectra 
in the individual J, H and K bands shows the continuum 
level increasing to longer wavelengths showing the develop- 
ment of an infrared excess associated with dust. 

The most interesting spectral features in V2615 Oph arc 
the first overtone CO bands in the K band which is discussed 
in coming sections. The other atomic lines in V2615 Oph are 
similar to those seen in carbon-oxygen novae and which were 
discussed in our earlier work on V1280 Sco - hence they are 
not discussed further here. There is, however, an unidentified 
line at 2.088 /xm which is discussed below. 

3.3 An unidentified line at 2.0888 /ttm: possibly a 
Fell line excited by Lyman a fluorescence? 

We note the presence of an emission line at ^ 2.0890 fim. in 
the K band that is seen fairly prominently in the spectra 
from 27 April onwards. We propose that this is a Fell line 
and additionally investigate whether this line could be ex- 
cited by Lyman a (Lya) fluorescence. In the near-infrared, 
there are a few Fell lines seen in the spectra of novae, which 
are believed to be primarily excited by Lyman a and Ly- 
man continuum fluorescence. Among these are the so-called 
"one micron Fe II lines" at 0.9997, 1.0171, 1.0490, 1.0501, 
1.0863, and 1.1126 /xm seen in several novae (Rudy et al. 
1991, 2000). In addition, two other Fe II lines at 1.6872 and 
1.7414 lira in the H band, are also proposed to be pumped 
by the same mechanism (Banerjee et al. 2009). The H band 
lines are prominently detected in the 2006 outburst of re- 
current nova RS Oph (Banerjee et al. 2009), in the slow 
nova V2540 Ophiuchi (Rudy et al 2002) and possibly also 
in the recurrent nova CI Aquila (Lynch et al. 2004). These 
H band lines could be present in the spectra of other novae 
too, but have evaded detection because of blending with the 



Table 3. List of observed lines in the JHK spectra 
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Br 11 (1.6806 jim) and CI lines at 1.6890 and 1.7449 ^im 
which lie close by. However these lines are well resolved in 
RS Ophiuchi (Banerjee et al. 2009), especially during the 
later stages of its outburst when all line widths in RS Oph 
become small - due to deceleration of the shocked, emitting 
gas - and blending effects are thereby minimal. 
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Figure 3. The J band spectra of V2615 Oph on different days witfi the fiux normalized to unity at 1.25 fim. 



The excitation mechanism of the Fe II lines appears to 
be a three step process (Johansson & Jordan 1984; Bautista 
et al. 2004 and references therein, Banerjee et al. 2009). Fell 
is first excited from low-lying levels by Lya or Lyc to a high 
energy level (typically 11 to 13 eV above ground state) which 
decays, in the second step, to feed the upper level associated 
with the observed Fell line (the decay of this upper level, in 
the third step, leads to the line proper). The 2.0888 ^m line 
resuhs from the decay of the 3d^{2F)4:S c^F term at ~ 6.209 



eV above the ground state. It so happens that this term can 
be fed by not just one, but in fact by several high lying levels 
- each of these high lying levels capable of being pumped by 
Ly a photons. As illustrative examples, it is noted that pho- 
tons at 1213.738, 1214.067, 1216.239, 1216.272 and 1217.152 
A, which are reasonably coincident with the Ly a line center 
at 1215. 67lA, can excite transitions from the low lying (5D) 
4s a^D term in Fe II (at 1.076 eV above ground state) to the 
higher excited levels ((5D) 4p 4F, (5D) 4p 4D and (3P) 4p 



1.5 



1.0 



28 Mar 




2.0 



1.0 



I .... I .... I .... I .... I .... I .... I . 



31 Mar 





2.0 



1.5 



1.0 



. I .... I .... I .... I .... I .... I .... I . 



02 Apr 




. I .... I .... I .... I .... I .... I .... I . 



2.0 



1.5 



1.0 



03 Apr 




05 Apr 



06 Apr 




1.5 



1.0 




V 








2.0 



> 



CD 



1.0 



3.0 



2.0 



1.0 



3.0 



2.0 



1.0 



07 Apr 



16 Apr 






3.0 - 



2.0 



18 Apr 







1.0 









27 Apr 




4.0 



2.0 - 



30 Apr 



02 May 




4.0 



2.0 




I .... I .... I .... I .... I .... I .... I . 



. I .... I .... I .... I .... I .... I .... I . 



. I .... I .... I .... I .... I .... I .... I . 



06 May 



08 May 



09 Jun 





2.0 



1.0 




2.0 



1.0 












I I 



1.50 1.55 1.60 1.65 1.70 1.75 1.80 1.50 1.55 1.60 1.65 1.70 1.75 1.80 1.50 1.55 1.60 1.65 1.70 1.75 1.8 

Wavelength (micron) 

Figure 4. The H band spectra of V2615 Oph on different days with the flux normalized to unity at 1.65 ^m. 



4P respectively at around 11.3 eV ). Since HI lines in no- 
vae (Ly ct included) are broad with widths extending upto 
a few thousands of km/s (1 Angstrom corresponds to about 
250 km/s at the Ly a wavelength), these photons could con- 
tribute to the flourescence process, even though they are not 
coincident with the Ly a line centre. These higher levels at 
11.3 eV can then decay via ultraviolet photons to the upper 
level of the 2.0888 ^m transition. The Kurucz atomic data 



Q on which the present analysis is based, indicates there are 
additional Lya fluorescing candidate lines (apart from the 
five discussed here), all within a few Angstroms of the Ly 
a line center - that could also contribute to the Ly a flu- 
orescence process. Therefore, the 2.0888 /xm line could be 



^ http:/ /cfa-www. harvard.edu/amp/ampdata/kurucz23/sekur. html 
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Figure 5. The K band spectra of V2615 Oph on different days with the flux normalized to unity at 2.2 /im. 



excited by Ly alpha flourescence, if its identification with Fe 
II is correct. 

A few cautionary words regarding the identification of 
the 2.0888 fj,m with Fell may be in order. In case the Fe 
II identification is correct, then a few other Fe II lines - 
as mentioned earlier - could also be expected viz. lines at 
0.9997, 1.0171, 1.0490, 1.0501, 1.0863, 1.1126, 1.6872 and 
1.7414 ^m. Unfortunately the first four of these lines are 
not covered in our spectra while it is difficult to make any 



definitive conclusion about the 1.0863 fim line which is in a 
region of low signal since it is at the edge of our spectral win- 
dow. However the 1.1126 fim line is seen. It is also difficult 
to draw a firm conclusion whether the 1.6872 and 1.7414 
^m lines are present. Unfortunately both these lines occur 
at positions close to strong CI and HI lines (Table 3) which 
could lead to line blending. In essence, further detections of 
the 2.0888 fim line in other novae is desirable to enable a 
secure identification. 
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3.4 Modeling the CO emission 

The CO emission is modeled by assuming tlie gas to be in 
thermal equilibrium characterized by the same rotational 
and vibrational temperature. The populations of the dif- 
ferent levels can then be determined from the Boltzmann 
distribution since the energy of individual rovibrational lev- 
els is known (we adopt a similar model as Spyromilio et al. 
1988). In our calculations, we have used values for the ro- 
tational and vibrational constants of ^'^C^^O, ^^C^^O from 
the NIST database Q except for the vibrational constants 
for 13 CO which are adopted from Benedict et al. (1962). 
We have not considered other isotopic species like ^^C^^^O, 
"C^'^O etc. since they are expected to have low abundances. 

The line luminosity E of each rovibrational transition 
can be calculated from knowing the population of the up- 
per level involved in the transition, the associated transition 
probability for the line (Goorvitch 1994) and the photon en- 
ergy (hu) associated with the transition. To enable a com- 
parison of the model value with the observed data, the line 
luminosity E calculated above, which is in units of ergs/sec, 
needs to be converted into units of the observed flux (we 
use units of ergs/s/cm^//j,m). This conversion is achieved 
by first dividing E by AtvD^ (where D is the distance to 
the source) and subsequently scaling i? to a unit of wave- 
length. In implementing the last step, it is assumed that 
each rotational hne is gaussian in shape. Then, the peak in- 
tensity of such a gaussian (plotted with its ordinate units in 
ergs/s/cm^//im; abscissa in units of /im) can be determined 
analytically by ensuring that the integrated area under this 
gaussian matches the known quantity E/AtvD^. The gaus- 
sian line profiles of all the rotational lines are thus generated 
and co-added together to generate the resulting envelope of 
the CO emission. To this envelope we add the appropriate 
continuum, the level of which is determined from broad band 
photometry of Table 2, to provide a model CO emission spec- 
trum for a particular day. Such model spectra matching the 
observed data of 28 and 31 March; and 2 and 3 April are 
shown in Figure 5. The parameters that are fed as inputs for 
the modeling are the total mass of the CO gas (A/co), the 
i2q^i3q j-a^io (designated as a) and the gas temperature 
(Tco)- Once Mco and a are chosen, the total number of 
i^CO and ^^CO molecules are fixed - thus the level popula- 
tions in thermal equilibrium at a temperature Tco can be 
calculated and the emerging spectrum determined. The es- 
timated model CO flux is therefore an absolute quantity but 
it could however be subject to certain errors. One of these 
is the value of the distance D to the source - the model CO 
flux scales as . 

Different combinations of the parameters Mco, ol and 
Too were tried to find the best fit between model spectra 
and the observed data. The characteristics of the CO emis- 
sion remains fairly constant between 28 March and 3 April. 
The profiles of Figure 6 are reasonably well modeled with 
closely similar values of Tco and Mco in the range of 4000 
- 4300K and 2.75x10"** to 3.25x10"* Mq respectively. Some 
of the difficulties encountered in the modeling, and also pos- 
sible sources of errors involved therein, are as follows. It is 
noted that changes brought about in the spectra by either 



^ http:/ /physics. nist.gov/PhysRefData/MolSpec 



changing Tco or changing Mco can be similar in the fol- 
lowing respect. Increasing either of these quantities enhances 
the absolute level of the CO emission. However, the contri- 
butions of these two parameters can be distinguished by the 
fact that increasing Mco just scales up the overall level of 
the CO emission level. On the other hand, increasing Too 
not only increases the CO flux but additionally changes the 
intensities of the different vibrational bands, relative to each 
other, within the first overtone. Therefore, vibrational bands 
from v^5, would help the analysis, but they are located in 
regions of poor atmospheric transmission, beyond the spec- 
tral coverage and the v=5-3 band is barely covered. 

While making the fits we note that the CO bands are 
likely contaminated by CI lines at 2.2906 and 2.2310 /xm 
and from Nal lines at 2.3348 and 2.3379 /im. The position 
of these lines are marked in the bottom panel of Figure 6 and 
there are discernible structures at these positions, in most 
of the profiles, indicating these lines are present. The 2.29 
/xm (v = 2-0) band would appear to be affected, especially 
by the CI 2.2906 line which can be significant in strength in 
carbon-oxygen novae (Rudy et al. 2003; Das et al. 2008) and 
whose presence is possibly responsible for lack of agreement 
between model and observed data in this region. In view 
of the above, we have relied more on the fits to the higher 
bands ( the u = 3-1 and 4-2 bands) while estimating the CO 
parameters. We have also assumed the CO emission to be 
optically thin. For the assumptions and constraints outlined 
above, the formal fits of Figure 6 yield estimates of Tco of 
4100, 4000, 4100 and 4300K (with an error of ± 400K) for 
28 March, 31 March, 2 April and 3 April. We find that a 
constant CO mass of 3.0x10"* Mq, or at most a marginal 
variation of Mco between 2.75x10"* to 3.25x10"* Mq, can 
account for the observed profiles on the four days. The fits 
shown in all panels of Figure 6 are made for the same mass 
of 3.0x10"* Mq and the implications of a fairly constant 
mass is discussed shortly. 

As in earlier studies on the ^^C/^^C ratio in a few no- 
vae, discussed in the following sub-section, we are also able 
to place only a lower limit on this ratio viz. ^'^C/^'^C is 
greater than 2. If this ratio is decreased below 2, the I'^CO 
contribution becomes rather prominent i.e. the ^^^CO band- 
heads begin to prominently appear in the synthetic spec- 
tra resulting in poor model fits. The determined lower limit 
for the i^C/^C ratio may be compared with that expected 
from theoretical calculations. For carbon-oxygen novae like 
V2615 Oph, the expected ^'^C/^'^C ratio in the ejecta has 
been computed (e.g. Jose and Hernanz (1998); Starrfield et 
al. (1997)) and shown to be dependent on the white dwarf 
mass among other parameters. Different models by Jose and 
Hernanz (1998) show that this ratio is approximately con- 
strained between 0.3 to 0.65 for a white dwarf mass between 
0.8 to 1.15 Mq; Starrfield et al. (1997) find the ^^C/^^C ra- 
tio to decrease from 2.4 to 0.84 as the white dwarf mass in- 
creases form 0.6 to 1.25 Mq. However, our observations and 
modeling seem to indicate that ^^C is possibly not synthe- 
sized to the high levels predicted by the theoretical models. 

To calculate the CO column density, we assume that 
the CO is uniformly spread in a shell of thickness AR 
and volume 47r7?^AR. The radius R of the shell is esti- 
mated kinematically knowing the time elapsed since out- 
burst and the velocity of the shell. A value of ~ 1450 km/s 
is adopted for the shell velocity based on the P-Cygni ter- 
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minal velocity reported by Munari et al. (2007) and simi- 
lar values of the FWHM of the near-IR lines as found in 
this work. The column density, which is proportional to 
(McoAR)/(47r7?^AR), is then independent of the thickness 
of the shell. In this manner we determine the CO column 
densities to be 7.5xl0^^ 4.1xl0^^ 3.8x10^® and 2.9x10^** 
cm~'^ on 28 March, 31 March, 2 April and 3 April respec- 
tively. 

It maybe useful to present similar results in other known 
novae in which the first overtone of CO has been detected. 
These are V2274 Cyg (Rudy et al. 2003), NQ Vul (Ferland 
et al. 1979), V842 Cen ( Hyland & Mcgregor 1989, Wich- 
mann et al. 1990, 1991), V705 Cas (Evans et al. 1996) and 
V1419 Aql (Lynch et al. 1995). The CO mass Mco was 
determined to be lo^s sio-s g 1.6x10"* Mq) in NQ 

Vul by Ferland et al. (1979) from observations 19 days after 
outburst or 20 days before the large visual fading associ- 
ated with dust formation. A limit of ^^C/^^C greater than 
3.0 and a temperature Tco = 3500 ± 750K were determined 
for the object. In V2274 Cyg, Rudy et al. (2003) from obser- 
vations 17 days after discovery, determined Mco ~ 8x10"^ 
Mq, i^C/i^C ^ 0.83 ± 0.3 and Tco = 2500K. In nova V705 
Cas, observations by Evans et al. (1996) taken at two epochs 
i.e. one day before maximum light and 26.5 days after max- 
imum light yielded estimates for Mco of (2.8 ± 0.2)xl0~^° 
Mq and (3.8 ± 0.2)xl0"^° M© respectively and Tco equal 
to 4300 ± 300K & 4500 ± 300K respectively. They estimate 
the ^^C/^'^C ratio to be J5 5 and the CO column density 
to be 2x10" cm"^ for Mco equal to 1.0x10"" Mq. Apart 
from the above sources, the ^^C/^^C ratio has also been esti- 
mated in two other novae. In DQ Her this ratio was found to 
be ^ 1.5 by Sneden and Lambert (1975); in V842 Cen obser- 
vations between days 29 and 45 after outburst by Wichmann 
et al. (1991) show the "C/"C ratio to be 2.9 ± 0.4. 

3.5 Evolution of the CO emission 

The detection of CO over a significant duration of time in 
the present observations, presents an opportunity - not avail- 
able before - to study the formation and destruction of CO 
during a nova outburst. The earliest theoretical studies of 
the chemistry of novae were done by Rawlings (1988) in the 
form of pseudo-equilibrium chemical models of the pre-dust- 
formation epoch. These models, which were developed with 
the main aim of explaining the observed presence of CO in 
novae, found that the outer parts of the ejecta have to be 
substantially more dense and less ionized than the bulk of 
the wind for substantial molecule formation to occur. For 
this to occur carbon has to be neutral. In a neutral carbon 
region, the carbon ionization continuum, which extends to 
less than 1102 A, shields several molecular species against 
the dissociative UV fiux from the central star. A more re- 
fined model for molecule formation in the nova outflow in the 
early stages is presented in Pontefract and Rawlings (2004; 
hereafter PR) - we try to correlate the present observational 
findings with the results in this work. The PR studies are a 
major extension of their earlier models with only one major 
qualitative point of departure viz. neutral-neutral reactions 
are now shown to be more important than photoreactions 
in determining the nova chemistry. 

A significant result in PR is the prediction of the evolu- 
tion of the fractional CO abundance with time. Two models 
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Figure 6. Model fits (dashed lines) to the observed first overtone 
CO bands in V2615 Oph. The fits are made for a constant CO 
mass of 3.0x10"® Mq on all days while the temperature of the 
gas Tco is estimated to be 4100, 4000, 4100 and 4300K (with an 
error of ± 400K) for 28 & 31 March and 2 & 3 April respectively. 
The bottom panel shows the position of certain CI and Nal lines 
that complicate the modeling. The other prominent lines seen in 
the spectra are Br 7 at 2.1655 /^m; NaT 2.2056, 2.2084 /im blended 
with weaker emission from CI 2.2156, 2.2167 /^m lines and other 
CI Unes between 2.1 and 2.14 ^m. The position of the i^cIBq 
bajidheads are .shown in the ton nanel. 
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are considered - model A considers oxygon rich ejecta and 
model B considers carbon rich ejecta. Figures 1 and 2 of PR 
show the evolution of the fractional abundance of different 
molecules and radicals, including CO, with time. It is seen 
that in both models, the CO abundance remains constant 
upto about 2 weeks after outburst (~ 12 days in case A and 
~ 15 days in case B). This behaviour, and the length of its 
duration, seems to be generic to the models. During this 
phase the CO is saturated - that is to say, all the available 
oxygen or carbon, whichever has the lower abundance, is 
completely incorporated into forming CO. After this there 
is a sharp decline in the CO abundance as CO is destroyed 
mostly by reactions with N and N+. During this stage, from 
Figures 1 and 2 of PR, we estimate an approximate decrease 
in CO by a factor of 1000 in ~ 27days for model A and a 
decrease by a factor of 100 in ^ 16 days for model B. The 
present data allows to check on these two vital aspects of 
the model predictions viz. the existence of a short-lived sat- 
urated phase followed by a phase involving rapid destruction 
of CO. 

The present observations and modeling show that the 
CO mass was constant between 28 March and 3 April i.e. for 
a period of 7 days. This puts a lower limit on the duration 
of the saturated phase since our observations commenced 
on 28 March, nearly 8-9 days after the beginning of the 
outburst on Mar. 19.812 UT. It is possible that the CO 
emission was present, and at similar levels, between March 
19 and March 28 also. After all, CO has been seen very early 
after commencement of the eruption as in the case of V705 
Cas (Evans et al. 1996). If that is the case, then the upper 
limit on the saturated-phase timescale would be around 15 
days. Thus the evidence indicates that a phase does exist 
when the CO mass is constant and whose duration ts is 
constrained within 7 ^ ts ^ 16 days. This observational 
finding conforms well with the predicted timescales of ~ 12- 
15 days from the PR model calculations 

Between 3 April and 5 April there is a sudden decrease 
in the CO strength. It is difficult to meaningfully model the 
5 April data, and hence derive the CO parameters, because 
the relative contribution of the CI 2.2906 /itm line is con- 
siderable at this stage. Also the signal-to-noisc ratio in the 
region around the CO emission region begins to get low now. 
However, if we assume that Tco has not changed drastically 
between 3 and 5 April, then from the diminished CO flux on 
the latter date, we estimate that Mco has decreased by a 
factor of 3. Beyond this date, the decrease in the CO emis- 
sion continues to take place, with a possible restrengthening 
again on 16 April, to finally drop below measurable limits by 
27-30 April. But the data beyond 5 April is of inadequate 
quality to make a quantitative assessment of the strength 
of the CO emission beyond this stage (or indeed whether 
CO is even present on some of the days; the presence of 
the CI lines complicates the assessment further). However, 
a quantitative assessment can certainly be made with a good 
degree of surety - there is a phase of rapid reduction in the 
CO emission after 3 April. The initial decrease is indeed 
rapid - the spectra over a 4 day gap, between 3 April and 
7 April, have only to be compared to establish this. This 
quantitative behaviour i.e. the rapid destruction of CO fol- 
lowing the saturated phase is again largely consistent with 
the predictions of the theoretical model (PR; 2004) . 

It is also possible to estimate the CO:C ratio, though 



with substantial uncertainties, at different stages of the CO 
evolution and compare it with model predictions. We as- 
sume the entire ejecta mass to be in the range of 10~^ to 
10~^ Mq which is fairly representative for the mass of novae 
ejecta. Carbon can be assumed to comprise about 10 percent 
of this mass as indicated by model calculations for elemen- 
tal abundances (Jose and Hernanz, 1998) in CO novae like 
V2615 Oph. Thus, in the initial saturated phase when Mco 
is found to be 3xl0~^ Mq, the CO:C ratio is determined to 
lie in the range between 10^^ and 10^^. This would clearly 
rule out the PR models in which the initial abundances of 
carbon are less than oxygen (model A) and be consistent 
with model B (carbon rich) where a CO:C ratio of ~ 10~^ is 
expected. At later stages, after the destruction of CO, Mco 
is difficult to estimate precisely but it is certainly lower than 
one-tenth of its value during saturation. This suggests that 
CO:C is in the range 10~^ to 10~^ and likely to be even 
lower. This is again reasonably consistent with the PR re- 
sults which find that in most models, regardless of whether 
the ejecta is O or C rich the CO:C ratio decreases with time 
and is less than 10~^ at a time greater than 50 days. 



4 SUMMARY 

Near-infrared spectroscopy and photometry of the dust 
forming nova V2615 Oph are presented. The key observa- 
tional result is the detection of first overtone CO emission 
in this nova. Modeling of the data indicates an initial phase 
when the CO is saturated and whose duration ts is con- 
strained in the range 7 ^ ts ^ 16 days. During this phase, 
the gas temperature and mass are found to be fairly constant 
in the range of 4000 - 4300K and 2.75x10"* to 3.25x10"** 
Mq respectively. A ratio of ^^C/^^C ^ 2 is inferred. Follow- 
ing the saturated phase, the CO is found to be destroyed 
fairly fast. The observed timescales involved in the evolu- 
tion of the CO emission and the estimated CO:C ratio are 
compared with model predictions and found to be in good 
agreement. 
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